format 8 

bind 
lambda   1.0
end

markov 10proc(c, mu) 
1 2 mu
1 0 lambda
2 1 2*lambda*c
2 3 mu
2 0 2*lambda*(1-c)
3 2 3*lambda*c
3 4 mu
3 0 3*lambda*(1-c)
4 3 4*lambda*c
4 5 mu
4 0 4*lambda*(1-c)
5 4 5*lambda*c
5 6 mu
5 0 5*lambda*(1-c)
6 5 6*lambda*c
6 7 mu
6 0 6*lambda*(1-c)
7 6 7*lambda*c
7 8 mu
7 0 7*lambda*(1-c)
8 7 8*lambda*c
8 9 mu
8 0 8*lambda*(1-c)
9 8 9*lambda*c
9 10 mu
9 0 9*lambda*(1-c)
10 9 10*lambda*c
10 0 10*lambda*(1-c)
end
10 1.0
end

* Functions & variables part:
func mtif(mu) mean(10proc;c,mu) / mean(10proc;c,0)


bind c 0.9999
expr c
expr 1.0/lambda, mtif(1.0)
expr 5.0/lambda, mtif(5.0)
expr 10.0/lambda, mtif(10.0)

bind c 1.0
expr c
expr 1.0/lambda, mtif(1.0)
expr 5.0/lambda, mtif(5.0)
expr 10.0/lambda, mtif(10.0)

end


